A one-level additive Schwarz preconditioner for 
a discontinuous Petrov-Galerkin method 



Andrew T. Barker 1 , Susanne C. Brenner 1 , Eun-Hee Park 2 , and Li-Yeng Sung 



1 A discontinuous Petrov-Galerkin method for a model Poisson 
problem 

Discontinuous Petrov-Galerkin (DPG) methods are new discontinuous Galerkin 
methods [3, 4, 5, 6, 7, 8] with interesting properties. In this article we consider a 
domain decomposition preconditioner for a DPG method for the Poisson problem. 

Let £2 be a polyhedral domain in W 1 (d = 2,3), be a simplicial triangulation 
of Q. Following the notation in [8], the model Poisson problem (in an ultraweak 
formulation) is to find <?/ G U such that 

b(^,f) = l(r) VreV, 
where [L 2 (Q)] d x L 2 {Q) xH?{dQ h ) xH~? {dQ h ), V = H(div;Q h ) xH l (Q h ), 
fo(*,y)= / oxdx — / uA\vxdx+ ^ / Cii-nds 



Y_, / c •gradvdx+ / vd„ds 

"en, JK Ken, JdK 



for % = (a,u,u, a„) e U and r = (t, v) e V, and / ( r) = J n fvdx. 

HereHQ /2 (dQ h ) (resp. H- l / 2 {dQ. h )) is the subspace of UKen h H 1/2 (dK) (resp. 
I\K£n h H~ l / 2 (dK)) consisting of the traces of functions in Hq(Q) (resp. traces 
of the normal components of vector fields in H(div;Q)), and //(div; £1/,) (resp. 
H l {£!,,)) is the space of piecewise //(div) vector fields (resp. H l functions). The 
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inner product on V is given by 

((ti,vi),(t 2 ,v 2 )) v = £ / [Ti -T 2 + divTidivT 2 + viv 2 +gradvi -gradv 2 ]c/x. 

The DPG method for the Poisson problem computes % G Uf, such that 

VreV,, (1) 

Here the trial space [//, (c I/) is defined by 

u h = n i^rfx n Pm(K)xp m+1 (dQ h )xp m (dQ h ), 

Ken h Ke£i h 
P„, (K) is the space of polynomials of total degree < m on an element K, P m+ i(dQ/,) = 

1/2 

H (oQf,) nUgsfii, Pm+\{dK), where P m+ \(dK) is the restriction of P m +\{K) to 
dK, and P m (di2,,) = H- X l 2 (d£2 h ) C\Y\ Kenh P m {dK), where P m (5/T) is the space of 
piecewise polynomials on the faces of K with total degree < m. 

Let V = {(T,v) G V : %\ K G [P m+2 {K)] d ,v\ K € P r (tf) VAT G X2 A } for some r > 
m + d. The discrete trial-to-test map 7), : £//, — !> V is defined by 

(r A *fc, r) v = b(vh,-r), V*;, G U h , r G V"', 

and the test space V/, is T^Uh. 

We can rewrite (1) as a/,^/,,, r) = l(Tf,r) for all w G t//,, where 

a h (?/,w) = b h {"//,T h v) = (T h *7/J h w) v 

is an SPD bilinear form on V/, x V/„ and we define an operator A/, : [//, — > U[ by 

(A/,*, r) = a A (*,*') V*,r G C//,. (2) 

Our goal is to develop a one-level additive Schwarz preconditioner for A/, (cf. [9]). 

To avoid the proliferation of constants, we will use the notation A < B (or B > A) 
to represent the inequality A < (constant) x B, where the positive constant only 
depends on the shape regularity of Qf, and the polynomial degrees m and r. The 
notation A « B is equivalent to A < B and B < A. 

A fundamental result in [8] is the equivalence 

ak(*,*) « llcrll^^ + IIbII^a) + II"IIhi/2 ( ^) + \\ 6 n\\ z H -V2 {dQh) (3) 

that holds for all % = (o,u,fi,O n ) G U/„ where 

II"IIhv 2(5%) = L II«IIhV2 W = L .if. .IMl£i W > w 

; Ken,, 1 KeQ h weH ( K )MdK=& 

\\ &tt \\ 2 H-W(dQ,) = £ W^l-ir-ldK) = E „,„ i nf , - ll«lll(div,*V (5) 
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Therefore the analysis of domain decomposition preconditioners for A/, requires a 
better understanding of the norms || ■ \\ H i/2/gg^ and || • \\ H -i/2^g K ^ on the discrete 
spaces P m+ \ (dK) and P m (dK). 



2 Explicit Expressions for the Norms on P m+ i(dK) and P m (dK) 

Lemma 1. We have 

IICII^ ) ^(llClli 2 (^ ) + I ICI^ (f) ) v;eP m+1 (dK), 

FEZ K 

where h K is the diameter ofK and Ek is the set of the faces ofK. 

Proof. Let JY{K) be the set of nodal points of the P m Lagrange finite element asso- 
ciated with K and jV(dK) be the set of points in JV(K) that are on dK. 
Given any £ G P m+ \{dK), we define g P m +\ (K) by 

?l . jUp) ifpe^{dK), 
Q*\p) = \ - (6) 
l£* ifpe^)\^(^), 

where is the mean value of £ over <9A". Since = £ on <9A", we have 

Suppose w £ H 1 (K) satisfies w = t, on <9/if . It follows from (6) and the trace 
theorem with scaling that 

l&lliw ^h K \\C\\l (d K) = MIHll 2 (<«) < \\w\\ 2 Hl{K) , (8) 
and, by standard estimates, 

IC*lirl(Jt) = IC* _ CdKlffl^ £ fyc K* ~ CdK\\l 2 {dK) 

= h- l \\w-w dK \\l i{dK) <\w\ 2 Hl{K) . (9) 

Combining (7)-(9), we have ~ \\t*\\ Z H \( K y The lemma then follows 

from (6), the equivalence of norms on finite dimensional spaces and scaling. □ 

Lemma 2. We have 
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Proof. We begin with the reference simplex K. Let RT m (K) be the ra-th order 
Raviart-Thomas space (cf. [2]). Given any £ £ P m (dK), we introduce a (nonempty) 
subspace RT m (K, \ ) = {q G RT m (K) : q ■ n = £ on dK and div q G P (K) } of RT m (K) . 
Let £, € ^^m^, C) be defined by 

qeRT m (K,Q) 2y ' 

Then the map S : P m (dK) — > RT,„(K) that maps £ to is linear and one-to-one, 
and we have (SQ • n = £ on div (SQ G P (K) and 

ll^llz,^) » IKIL 2( ^) y^UdK). (io) 

Let £1,...,^ be a basis of P,„(dK) and 1 = 0i,...,0jv„, G H l l 2 (dK) satisfy 

det [ £0,- ds] yt 0. We define the map Q : //(div; K) — > P,„(dK) by 

L J l<i.j<N,„ 

(Qq)$jds = (q ■n,0j) ff -i/2 (aje)xH i/2 (3J? - ) for 1 < ; < JV m . 

./ <7/l 

It follows from the definition of Q that ||2<7|Il 2 (,9,k~) S I Ml //(div it) f° r a ^ 1 e 
//(div; K), and £2<7 = £ if q ■ n = £ G P m (dK), in which case 

\\S£\\l 2 (K) ~ IKIIi 2 (5if) = ll&llz^*) £ ll?llff(cfiv;*)- C 11 ) 

Moreover, since 0i = 1, we have 

ldiv(SQdx = (Qq)lds=(q-n,l) H -i/2 {g&)xH i/2, d g :) = Idivqdx 

J K J oK J K 

and hence 

l|div(5C)|| L2 ^<!|div ? || i2 ^. (12) 

Now we turn to a general simplex K. It follows from (10)— (12) and standard 
properties of the Piola transform for //(div) (cf. [10]) that there exists a linear map 

5 : P m (dK) — > RT,„(K) with the following properties: 

(i) (SQ • « = C an d hence 

KWH-WdK) = , i° f Jtfllff(div;Jr) < \\SQ H (6iv;K) V£ € P m (dK), 

(ii) for any g G //(div; /T) such that q ■ n = £ , we have 

1 1 SCI I //(div; tf) < ||?||#(div;is:); 

(iii) div (5C) G P (K) and hence 

[diw(SQdx=[ Cds or ||div(5C)||l 2W = ( / C^)Vl, 
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(iv) we have 

/^llscilL^ W^iicilW 

Properties (i)-(iv) then imply 

K\\W [dK) « ll s CII/W) ^h K K\\l id K)+hl d (l dK Cds) 2 . □ 



3 A Domain Decomposition Preconditioner 

Let Q be partitioned into overlapping subdomains Q\ , . . . , £2j that are aligned with 
Q.fi- The overlap among the subdomains is measured by 8 and we assume (cf. [11]) 
there is a partition of unity 9\,...,9j G C°°(Q) that satisfies the usual properties: 
Oj > 0, Ej=i 0j = l on Q, 9j = on Q \ % and 

l|V0y|| Mfl ) < 5- 1 VI <j<J. (13) 

We take the subdomain space to be Uj = { >// £ U/, : * = on Q \ Qj}. Let w = 
(<J,u,G, d n ) G i//,. Then G t/ ; - if and only if (i) a and u vanish on every K outside 
£2j and (ii) u and a„ vanish on dK for every K outside Qj. We define «;(•,-) to be 
the restriction of «/,(-, •) on Uj x f/j. Let A 7 : Uj — > Uj be defined by 

(Aj^,wj) = aj{%), Wj) y^,wj G Uj. (14) 

It follows from (3) that 

aji-wj^j) » ll^-lli^.) + H^ll^^.) + ll"Jtfi/2 ((9%A) + ll^jllff-iA^)' ( 15 ) 

where = (<Jj,Uj,Uj, d n j) G f/y, fl.-fl is the triangulation of Qj induced by Q/, and 
the norms || • ||#i/2oq. ) and || • \\g-i/2ig a . ) are analogous to those in (4) and (5). 

Let /y : {7; — > Uh be the natural injection. The one-level additive Schwarz pre- 
conditioner Bjj : U' b — > Uh is defined by 

B k = ilAj ll 'r 

y=i 

Lemma 3. We have 

Ki n {B,,A h )>8 2 . 

Proof. Let If, i, //, 2, ^,3 and //, 4 be the nodal interpolation operators for the compo- 
nents l\K££2 h [ p m(K)] d , UKen,, p m (K), P m +\ (dil h ) and P m (dQ h ) of U h respectively. 
Given any % = [a,u,u, a n ) G {//„ we define %j G Uj by 

*j = {h,i(0jO),I ht 2(9ju),I h> 3{Gja),I ht4 (9j& n )). 
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Then we have % — Yjj=\ *? an d, in view of (14) and (15), 

(Aj^) a 114,1(0^)11^^.) + ||4. 2 (0 7 -«) ||i 2(fl . } 

+ \M°J*)\\Wp aj j + \M e J^)f B -l H9a ^y (16) 

The following bounds for the first two terms on the right-hand side of (16) are 
straightforward: 

IIM^)llL(^)<!l< 2 ( % ) and 114,2(^)11^.) <\Hl 2{£lj y (17) 
We will use Lemma 1 and Lemma 2 to derive the following bounds 

llA^^II^/a^^s- 2 !!^!^^, (18) 
II^^^^OII^,/,^,^^ 2 !!^!^^^^,. (19) 

Let K £ Qj.h- It follows from Lemma 1,(13) and standard discrete estimates that 

l|4,3(0^)|| 2 ? i/ 2w «fe(||4,3(0^)||i 2 ( 3jr) + E MOjQ&ip)) 

FeE K 

<hKMl {dK] +hK L (IIVflyllL^llflll^ + ll^llLwKifF)) 

^^ll"lli 2 0^)+ / ^ 5_2 ||"lli 2 (^)+ / ^ L l fl l|i(jn ~ 5_2 ll"llfli/2( 3 ]c)- 

Summing up this estimate over all the simplexes in Qj /, yields (18). 
Similarly, it follows from Lemma 2 and (13) that 

||4,4(0A)||ff-i/2 (3Jf) ~h K \\I hA (e j 6 n )\\l 2{dK) +h- d ( J I hA (9jd„)ds} 

% hK\\On\\ 2 Ll(dK) + hKd- 2 \\a n \\ 2 L2{dK) + h~ d ( ^ d„ ds) 2 < 8- 2 \\d n f H _ 1/2 m , 

where 0y is the mean value of Gj over K. Summing up this estimate over all the 
simplexes in Qjh gives us (19). 

Putting (2), (3) and (16)-(19) together we find £j =1 (Ay^,«?) < 5 _2 (A A *,*), 
which implies A m i n (B/ J A/ ? ) > S 2 by the standard theory of additive Schwarz precon- 
ditioners [11]. □ 

Combining Lemma 3 with the standard estimate A max (fi/iA/j) ^ L we obtain the 
following theorem. 
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Theorem 1. We have 

(n A s _ K^(B h A h ) 2 

Anin(,o/i^/ij 

where the positive constant C depends only on the shape regularity of and the 
polynomial degrees m and r. 

Remark 1. Theorem 1 is also valid for DPG methods based on tensor product finite 
elements. 



4 Numerical results 

We solve the Poisson problem on the square (0, l) 2 with exact solution u = 
sin(7Dci) sin (jd^) and uniform square meshes. The trial space is based on Q\ poly- 
nomials for a and u, P2 polynomials for u, and Pi polynomials for a n . We use 
bicubic polynomials for the space V in the construction of the trial-to-test map 7},. 

The number of conjugate gradient iterations required to reduce the residual by 
10 10 are given in Table 1 for four overlapping subdomains. The linear growth of 
the number of iterations for the unpreconditioned system is consistent with the con- 
dition number estimate Jc(A/,) < h~ 2 in [8]. Note that in this case the boundary of 
every subdomain has a nonempty intersection with d£2 and it is not difficult to use a 
discrete Poincare inequality to show that the estimate in Theorem 1 can be improved 
to K(BhA/ 1 ) < \ \nh\S~ l . This is consistent with the observed growth of the number 
of iterations for the preconditioned system as 8 decreases. 



Table 1 Number of iterations for the Schwarz preconditioner with subdomain size H = 1/2. 



h S I unpreconditioned preconditioned 



2- 2 2- 2 


496 


14 


2 -3 2 -3 


1556 


17 






14 


2 -4 2 -4 


3865 


20 


2- 3 




17 


2- 2 




14 


2- 5 2- 5 


8793 


27 


2- 4 




20 


2- 3 




18 



In Table 2 we display the results for h = 2~ 5 and various subdomain sizes H with 
S = H/2. The estimate fc(B/ ; A/,) < S~ 2 w H~ 2 is consistent with the observed linear 
growth of the number of iterations for the preconditioned system as H decreases. 
Such a condition number estimate for the one-level additive Schwarz preconditioner 
is known to be sharp for standard finite element methods [ 1 ] . 



s 
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Table 2 Number of iterations with h = 2 5 and various subdomain sizes H with 5 = H/2. 



h H |unpreconditioned preconditioned 



2- 5 2- 1 


8793 


15 


2- 2 




25 


2- 3 




45 


2 -4 




89 
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